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Abstract. The zero-temperature XX chain is studied with emphasis on the 
properties of a block of L spins inside the chain. We investigate the quantum 
fluctuations resulting from the entanglement of the block with the rest of the chain 
I using analytical as well as numerical (density matrix renormalization group) methods. 

$Zj ' It is found that the rest of the chain acts as a thermal environment and an effective 

temperature can be introduced to describe the fluctuations. We show that the effective 
temperature description is robust in the sense that several independent definitions 
■ (through fluctuation dissipation theorem, comparing with a finite temperature system) 

^ , yield the same functional form in the limit of large block size {L — s- oo). The effective 

' temperature can also be shown to satisfy the basic requirements on how it changes 

. when two bodies of equal or unequal temperatures are brought into contact. 

a^ 
o 

^ . 

^ ■ 1. Introduction 

Entanglement properties of various spin chains have attracted much attention recently. 
. The most intensively investigated cases are spin chains in their zero temperature ground 
C . states, where the entanglement between a contiguous block of spins and the rest of the 
Q . chain is measured through the von Neumann entropy. Thus far, the entanglement 
entropy has been found to be a well behaving measure of entanglement which shows 
^ ! some universal features at quantum critical points [B IS |3l HI [5]. Despite capturing 
some important aspects, one should not forget, that we inevitably lose information on 
the entangled state by using only a coarse grained measure. Indeed, the restriction of 
the pure ground state to a subsystem results in a mixed state with quantum fluctuations, 
and the weights are given by the full spectrum of the reduced density matrix, which 
means a number of parameters proportional to the exponential of the system size. Hence 
an "integrated" measure such as the entropy will not tell about the distribution of the 
fluctuations. This gives an obvious motivation for a deeper analysis of reduced density 
matrix spectra. 
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In the past few years, much effort has been devoted to the investigation of reduced 
density matrices in a series of simple model systems O O [HI El [lO]- The results 
obtained indicate that, in several cases, the density matrix itself can be expressed as the 
exponential of a quadratic Hamiltonian. An example is the case of free fermions in one 
dimension (equivalent to an XX chain) which was investigated in detail pTl [T2l [13] and 
the reduced density matrix for arbitrary subsystem sizes has been obtained in a closed 
form. 

Here we try to develop a simple interpretation of the above results, concentrating 
on the XX chain. Looking at a subsystem of L spins in an infinitely large chain, the 
quantum fluctuations in the block are due to the coupling with the rest of the chain. 
One could think of an analogy with the subsystem-reservoir scheme of the canonical 
ensemble in statistical physics where now the spins outside the block would play the 
role of the reservoir. Therefore it is an intriguing question whether it is possible to define 
an "effective temperature" to describe the zero temperature quantum fluctuations with 
Boltzmann-like distributions, meaning that the weights of the reduced density matrix 
Pl decay exponentially with the subsystem energies: 



inverse temperature. The aim of this paper is to show that there are several ways to 
define an effective (inverse) temperature but they are all consistent in the sense that 
the leading order dependence on the subsystem size is identical, given by: 



where J is the coupling between nearest-neighbour spins. In the following we shall 
present a detailed derivation of the above result as outlined below. 

In Section [21 we give a brief summary of the results on the reduced density 
matrix of the XX chain |TT]. Next (Sec. [3]), we analyze the one-particle spectra of 
the effective Hamiltonian emerging in the formula of the density matrix and that of 
the Hamiltonian describing a finite XX chain. We define the effective temperature 
starting out from a continuum limit formula for the eigenvalues [13], then verify these 
results by analyzing the multiparticle spectra of finite systems either by means of the 
operator based truncation scheme [12] or density matrix renormalization group (DMRG) 
calculations [T4j . 

Section [H is concerned with an alternative definition of the effective temperature, 
using an extension of the fluctuation dissipation theorem to zero temperature. The 
calculations are carried through for two different conserved quantities and the results 
obtained show agreement with the previous definitions. 

With the concept of a system in an effectively thermodynamic environment, we try 
a direct mapping in Section [51 The effective temperature is determined numerically 
by comparing canonical and zero temperature expectation values. Surprisingly, we 




^^"^ = ^h^ 



(2) 
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recover the results of the fluctuation-dissipation ratio with a very high precision, and 
an explanation for the precise agreement is also given. 

The definitions are also tested by studying their consistency with thermodynamics 
(Sec. |6]). Namely, we ask what happens when subsystems of equal or different 
temperatures are brought into contact. We find that equal temperatures do not change, 
while in case of unequal temperatures the resulting temperature is in between the initial 
ones. 

Finally, we summarize our results (Sec. [T]) emphasizing also the problematic points 
which might set a limit on the applicability of the simple picture described here. 



2. Reduced density matrix for the XX chain 

Our starting point for the calculations is the Hamiltonian of the XX chain 

N-l N 

i=i i=i 
where s"(a = x, y, z) are the Pauli spin matrices at sites j = 1,2, . . . , N of the chain, J 
is the coupling and h is the magnetic field. All our results, unless stated otherwise, are 
derived in the case of zero magnetic field {h = 0). It is well known that this model can 
be transformed into a chain of free Fermions with Hamiltonian 

where and Cj are fermionic creation and annihilation operators and energy is measured 
in units of the coupling J. 

If one considers only a subsystem of size L in a large chain with spins, the state 
of the subsystem is fully described by its reduced density matrix = tr7v-L|^o) (^ol; 
obtained from the ground state |\E'o) by tracing out external degrees of freedom. It was 
previously shown [11] that for a chain with Hamiltonian (jl]) the reduced density matrix 
can be written in an exponential form 

h = ie-^ (5) 

where Z = Tr(e~'^) is a normalization constant and H is an effective Hamiltonian of 
the following form: 

H = J2A^Acr (6) 

One can see that H is quadratic in the fermionic operators. The matrix elements Aij 
can be given through the elements of the correlation matrix Cij = (cjcj) by the following 
relation: 

^ = ln^, (7) 
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where A and C denote the matrices composed of the elements Aij and Cjj, respectively. 
Hence, the effective Hamiltonian can be diagonalized, and the eigenvalues can be 
obtained through the eigenvalues ^ of the correlation matrix: 

H = Y.euflfk. e. = lni^, (8) 

where fl and fk are fermionic operators obtained by unitary transformation. 

The diagonalization of the correlation matrix is especially simple if we take the 
whole system to be infinitely large (A^ oo). In this case the elements of the matrix 
depend only on the distance r = |i — j| between the sites and can be written as: 





, r odd 






0, 


r even 


1 

2' 


r = 



sin(f(2-j)) _ ... 
L-^j - 7r(i-j) ~ > 0' reven • 

For finite N , each element of the correlation matrix has to be evaluated numerically as 
a sum of terms. 

The purpose of studying the structure of the reduced density matrix came from 
our intention to define an effective temperature by making the correspondence between 
the density matrix weights and Boltzmann factors ([1]) with respect to the subsystem 
Hamiltonian II i. This latter one can be diagonalized exactly and it also transforms to 
a system of noninteracting fermions: 

^L = ^Afc4cfc, Afc = - cos ( ^ j . (10) 
k=l V + / 

Peschel's results show that pl can be written as an exponential of a free fermion 
Hamiltonian ([8]). Thus, from comparison with (ITOl) it can be seen, that all the 
information we need is stored in the one-particle spectra Ek and Afc, where the former 
one has to be numerically calculated. 



3. One-particle and multiparticle spectra 

Figure [H shows the one-particle spectra of the effective Hamiltonian H for various 
subsystem sizes with the chain size N taken to be infinite. First of all, we should 
notice the apparent particle-hole symmetry of the spectra, which is made more visible 
by shifting the curves to a common origin. This symmetry property is not surprising 
as it is already present in the spectrum of the XX chain with zero magnetic field. The 
saturation of the curves for \k — {L + 1)/2| > 10 is due to numerical errors, since the 
major part of the eigenvalues of the correlation matrix Cjj lie exponentially close to 
and 1. 

The weights of the density matrix can be produced from these one-particle spectra in 
the same way as in the case of a real excitation spectrum: the largest weight is obtained 
by filling the modes up to the Fermi level, i.e. all modes that are negative. The next 
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Figure 1. One-particle spectra of the effective Hamiltonian H for various block 
sizes L with the chain size N taken to be infinite. The spectrum is shifted for 
each L to visualize the particle-hole symmetry. The saturation of the curves for 
|fc — (i + 1)/2| > 10 is due to numerical errors. 



largest weights are obtained by filling the smallest positive Sk ( "particle excitation" ) or 
removing the smallest negative one ("hole excitation"). Thus, the problem of finding 
the 2^ weights of the density matrix can be reduced to the diagonalization of an L x L 
matrix. 

3.1. Continuum limit 

Taking a continuum limit in the block {L oo) Peschel calculated analytically the 
lowest lying eigenvalues [T3] : 

One can see that the spectrum is linear near the Fermi level and the eigenvalues scale 
with 1/lnL. Looking at the excitation spectrum of the block XX Hamiltonian H^, 
one finds that it is also linear near the Fermi level and, for large block sizes, it can be 
approximated as: 

/ n n 2j — 1\ TT, , 



Fluctuations in the zero temperature XX chain: Effective temperature 



6 



Hence, the effective inverse temperature can be obtained by comparing only the low 
lying part of the one-particle spectra. From Eq. ffTTl) and f|T2|) one has: 

A« = (13) 

It can be seen that, according to the expectations, /3e// — >■ oo when L — > oo, that is we 
recover the zero temperature ground state if we take the whole chain as the subsystem. 
It has to be emphasized that the above definition of the effective inverse temperature 
describes correctly only the low energy part of the spectrum. However, it is well known 
from DMRG calculations that indeed, it is only a small portion of the weights which are 
relevant for the physics of the system |14j, therefore we also expect that the effective 
temperature is reasonably defined. 

We should note that, comparing with numerical results, one finds deviations from 
the asymptotic formula ( fT3i) . In the paper of Cheong and Henley [12], where the idea 
of defining an effective temperature has first appeared, the logarithmic scaling of the 
eigenvalues could not be detected for system sizes L < 100, hence they found that 
Peff oc L. This can be understood since the logarithmic dependence on the block 
size makes convergence very slow and one needs corrections even for large block sizes 
(L = 2000 — 5000). Determining Peff as the ratio of the lowest eigenvalues of ( fTTi) and 
(fT2|) we have tried a fit a x L/(lnL + b) which gave values a = 3.12 ^ it and b = 2.58 
which shows considerable finite size corrections, already noted in [13] . 

3.2. Operator-based truncation and DMRG results 

In order to compare the results of the continuum limit for the effective temperature 
( |T3l) to the one defined in Eq.([l]), one has to consider finite size systems and build 
the multiparticle spectra from the one-particle spectra. A possible method is the so- 
called operator-based truncation scheme which was devised by Cheong and Henley [12] 
and consists of truncating the one-particle spectrum by retaining only a small number 
of the excitations symmetrically around the Fermi level. The method can be applied 
successfully, since the normalization constant Z itself can be determined by using only 
the retained eigenvalues. Indeed, it can be verified [12] that 

1 ^ 1 

^ = det(l - C) = n . (14) 

k=l 

For large positive eigenvalues 1 + e~^'= ~ 1 while for large negative ones 1 + e~^* ~ e~^* 
thus the latter terms simply cancel the factors of the occupied one-particle states in the 
exponential ([5]). Hence the most relevant weights of the reduced density matrix can be 
reproduced from the low lying excitations of the one-particle spectrum and the highly 
excited states, corrupted by numerical error, can be discarded. 

The same method can be applied to obtain the multiparticle energy eigenvalues of 
the block, except that in this case we can calculate the given subset of the spectrum 
exactly, since the one-particle energies are known analytically. With this systematic 
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procedure a one-to-one correspondence between the block energies and the density 
matrix weights can be estabhshed. 

The multiparticle spectra constructed from the one-particle spectra has also been 
confirmed by the DMRG method [14j . In the DMRG method, the finite system is divided 
into two subsystems, thus the multiparticle spectra uJa^L) of the reduced subsystem 
density matrix and the eigenvalues Ea{L) of the reduced subsystem Hamiltonian can 
be calculated directly for a block of L spins, where a runs from one to the number 
of basis states, M, used to describe the block. We have calculated u!a{L) and Ea{L) 
as a function of the subsystem size L on finite chains with lattice sites with open 
boundary condition using the dynamic block-state selection (DBSS) approach [16]. The 
DBSS method allows a more rigorous control of the numerical accuracy by setting the 
threshold value of the quantum information loss x ci^nd the minimum number of basis 
states, Mmin, used to describe the subsystem blocks a priory of the calculations. For a 
given finite system with N lattice sites subsystem blocks with L = 2, . . . , N — 2 can be 
obtained using the sweeping procedure of the so called finite lattice algorithm. 

First, we have compared exact results obtained by exact diagonalization for L = 
N/2 on a short = 20 chain to the DMRG results by systematically adjusting x- In 
Fig.[2]we plotted uJa{L) as a function of Ea{L) obtained after the sixth DMRG sweep for 
various x values with M^^^ = 16. The accuracy of the diagonalization of the superblock 
Hamiltonian was set to 10~^. The agreement between the exact diagonalization and 
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Figure 2. Eigenvalue spectrum of pL as a function of the eigenvalues of Hl with 
i = 10 for a chain with = 20 lattice sites for various x values. The result of the 
exact diagonalization is shown by the solid line. 



DMRG results is evident from the figure for large uJa{L) even for large x values. On 
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the other hand, for the small uJa{L) spectrum more accurate calculations with smaller 
X values have to be performed for better agreement. 

In order to determine the size dependence of the effective temperature similar 
calculations were repeated first for subsystems with L = N/2 by varying N and then 
by keeping the size of the system fixed but varying the size of the block. For the 
former case Fig. [3] shows our results for = 50, 100, 200, 400, 800 using DMRG method 
with six sweeps, x = 10^*^, Mmin = 1024 and the operator-based truncation method 
with 20 particle-hole excitations kept. For large uJa{L) eigenvalues the agreement 
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Figure 3. Eigenvalue spectrum of pL as a function of the eigenvalues of Hl with 
L = N/2 for N = 50,100,200,400,800 using the operator-based truncation method 
and the DMRG algorithm with six sweeps. Degeneracies in the spectra are removed 
and all Ea{L) are measured from zero. 



between the operator-based truncation and DMRG result is evident from the figure, 
while for the lower part of the spectrum the deviation between the two approaches due 
to truncation procedures is more significant. This result, however, confirms the validity 
of the operator-based truncation method and it is sufficient in order to determine the 
effective temperature. Data points corresponding to uJa{L) > 10~^ eigenvalues for each 
L can be fitted with a straight line in order to obtain (3eff{L) which is shown in Fig. HI 
Our result indicates that (3eff{L) can be fitted in leading order with the following ansatz 
a X L/ {In L + b) and we obtained 6.42 > a > 5.76 and 2.12 > b > 1.57 depending on the 
number of data points taken into account when (3eff{L) was fitted by the straight lines. 
Using the operator-based truncation approach we repeated the same calculations with 
L = N/A for N = 200, 400, 800, 1600, 3200 and found again that a changed between 6.37 
and 5.63. The value of a close to 27r is twice that of given in Eq. (fT3!l . This is due to the 
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fact that the blocks were always determined for a semi-infinite system having one open 
boundary. When the same calculations are carried out for chains with up to = 200 
lattice sites using DMRG with periodic boundary condition (PBC) in which case there 
are two couplings across the block interface we obtained 3.3 > a > 2.9 and 2.1 > b > 1.8 
in agreement with Eq. (JT3i) . 

400 1 ' ' ' ' ' ' ' 1 

350 



ca 




50 100 150 200 250 300 350 400 

L 

Figure 4. f3eff{L) as function of i = N/2 obtained as the absolute value of the slopes 
of straight line fits to the ln{u!a{L)) vs. Ea{L) data sets. The solid line is the result 
of the fit with a = 6.309 and b = 2.101. 

We have also calculated uJa{L) and Ea{L) for various finite chains with = 
50, 100, 200, 400 lattice sites for DMRG blocks with L = 2, N/2 using several sweeps of 
the finite-lattice algorithm. For a given N taking data set corresponding to uJa{L) > 10~^ 
a linear fit for each L can be obtained in order to determine f3eff{L) as a function of 
L. Our result shown in the left panel of Figj5] indicates that in leading order f3eff{L) 
can again be fitted very well with the following ansatz a x L/(lnL + 6). For all fits we 
obtained 6.38 > a > 6.21 and h changed between 2.08 and 2.49. We also found that 
these fits are more stable than the above described previous approach. The value of a 
is again very close to 27r. When the same calculation is repeated using PBC we found 
3.28 > a > 2.94 and 2.12 > 6 > 1.81 in agreement with Eq. (fT3|) as also shown in the 
right panel of FigJSl The factor of two in a for OBC and PBC is the same effect that 
has been discussed for the entropy profile of a finite subsystem for critical models with 
open or periodic boundary conditions |4j. Finally, we conclude that a straight line fit 
is not enough to describe Peff{L) as a function of L, but logarithmic corrections still 
could not be confirmed rigorously due to the limited system sizes. 
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Figure 5. Peff{L) as function of L obtained as the absolute value of the slopes of 
straight hne fits to the ln{u}a{L)) vs. Ea{L) data sets corresponding to finite chains 
with TV = 50, 100, 200, 400 lattice sites. DMRG results, shown in the left and right 
panel, were obtained with open and periodic boundary condition, respectively. The 
solid lines are the results of the fitting. 



4. Fluctuation-dissipation ratio 

To support the idea of an effectively thermodynamic system, it would be desirable to 
find some alternative way of defining a temperature, and to compare the results of 
the different approaches. In ordinary thermodynamic systems, coupled to a reservoir 
at some finite temperature, the fluctuation-dissipation theorem gives a possible way of 
defining the temperature. If we restrict ourselves to the static limit and consider only 
conserved quantities of the Hamiltonian, it reads: 

XA = (3{{A'),-{A)l), (15) 

where A is a conserved quantity, and xa = -^f- is the response function of A to an 
infinitesimal conjugate field Ha- The brackets denote averages taken with the canonical 
ensemble. 

The straightforward generalization to our case would be to define the effective 
inverse temperature as the ratio of the response function to the fluctuations inside 
the block. We shall present here the calculations for two cases, using the transverse 
magnetization or the magnetization flux as the observable A. From now on we shall 
consider only the case when the chain size is taken to be infinite. 
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Fluctuations of the transverse magnetization 

The transverse magnetization is a conserved quantity in the XX model, thus it is a good 
candidate for the definition of the effective temperature, as described above: 

where xm is the magnetic susceptibihty. Ml is the transverse magnetization of the block, 
and the brackets now denote averages taken with pL. It is important to observe that 
we do not need the reduced density matrix to evaluate the expectation values since it 
is equivalent to calculating them in the ground state |\E'o) of the whole chain, which is 
known exactly. 

Let us begin by calculating the block magnetization in case of an applied magnetic 
field. The ground state \^h) is known exactly [15], and we have: 

{^k\Ml\^k) = X](c]c, -\) = ^r (4c,)dg - § = ^(f - (17) 

where qh = arccos(— /i) is the Fermi wave number. The magnetic susceptibility is 
therefore 

L 



d{ML) 

Xm = 



dh 



(1^ 

TT 

h=Q 



The fiuctuations in flTB]) have to be evaluated at /i = 0, when the block magnetization 
(Ml) is vanishing, hence the fiuctuations read: 

{Ml) = i^^sls]) = (4c.c]c,) - ^. (19) 

Using Wick's theorem and introducing the difference r = i — j oi the sites we obtain 

(AfE) = -2t = + (20) 

i<j=2 r=l 
r odd 

where the elements of the correlation matrix Cij were given in Eq. ([9]). It is easy to see 
that the linear terms in L cancel, and the asymptotic form of the fiuctuations will be: 

(M^) ^ i- (In L + c), c = ln2 + 7 + l. (21) 

Note, that (12T!) is consistent with a result in [17] where the leading order size dependence 
of the particle number fiuctuations has been obtained in case of free gapless fermions in 
arbitrary dimensions. Therefore the effective inverse temperature is obtained as: 

Peff = ^r^— , (22) 

m L + c 

which is a form very similar to (fT3l) . however there is a correction to the logarithm 
already seen at the numerical evaluation of the one-particle spectra. 
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4-2. Fluctuations of the magnetization current 

To further verify the vahdity of our description, we shall calculate the fluctuation- 
dissipation ratio for an other conserved quantity, the magnetization current. When the 
whole system is in equilibrium, the expectation value of the current flowing through 
the subsystem is vanishing, though there are still current fluctuations inside the block, 
which can be calculated using the equilibrium ground state. 

However, in order to obtain the response function, one needs to drive the system out 
of equilibrium. Therefore we need the ground state of the chain with a magnetization 
flux, a problem which has been solved by Antal et al in [18] . We briefly summarize the 
method leading to the main results and refer to [H] for the details. 

First, one needs to calculate the current density 

1 = sfsf+i - sfsf^i = ^(4q+i - 4_^iq), (23) 

and then introduce a modified Hamiltonian H^^ = H ~\J with the Lagrange-multiplier 
A, where J = 'YliJi- Note, that A will now play the role of the conjugate field of the 
magnetization current. 

The ground state of the above modified Hamiltonian will give us a current carrying 
steady state of the XX chain with a prescribed amount of current. Since [H^ J] = 0, 
Hamiltonian H^^ can be diagonalized exactly and the resulting one-particle spectrum 
has the following form {h = 0): 

A, = —cos{q-^), (24) 

COSip 

with (f = arctan(A). One can see that the spectrum is just the original XX spectrum 
shifted by the wavenumber ip, and rescaled by cosip. Having obtained the structure 
of the ground state I^E'a), the elements of the correlation matrix can be evaluated and 
yield: 





-e"'^, rodd 


Trr 




0, 


r even 


1 

2' 


r = 



J. I nr ' 

{^xlqcml'^x) = g{m~l) = i 0, reven , r = m-l. (25) 



We are interested in the response function of the block to an infinitesimal change 
in the conjugate field A. The expectation value of the magnetization current has to be 
evaluated in the current carrying ground state: 

{^x\Jl\^x) = ^J2 (^(1) - = ^ sin^ = ^"^ J— , (26) 

2Z TT TT a/1 + A2 

hence the response function is given by: 



aw 



(27) 



TT 
A=0 



which, for large values of L, is the same as the magnetic susceptibility 0181) . 
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Next, the current fluctuations have to be calculated in equilibrium. With {Jl) = 0, 
we need only: 

1 

01) = -4 XI (^4^1+1 - cJ+iQ)(4c™+i - cl^^Cm))- (28) 

l,m=l 

Using Wick's theorem, the contribution from terms with \m — l\ > 1 is 

E [\9{m-l)\'-Re{g{m-l + l)g*{m-l-l))], (29) 

|m-/|>l 

while the |m — / | = 0, 1 terms give 



L-1 



g(0){l-g{0)) + 



g\l)+g\-l) 



(L- 2)1^(1)11 (30) 



2 

Introducing r = m — I and using the equilibrium value of the matrix elements g{r), the 
fluctuations read: 

L-1 fl 1\ 1 f^L-l-r ^L-l-r\ , , 

^ ' \r even r odd / 

The asymptotics of the above expression can be given as 

~ 4(ln^ + c'), c' = ln2 + 7 + l/2, (32) 

which differs slightly from (pT!) in the next to leading order term. Looking at Eqs. (pTll 
and fl32|) we recover again the same leading order size dependence of /5e//. 

To conclude, we emphasize that the results obtained for the effective temperatures 
coincide in leading order. This is important since now contributions from all of the 
excited states are included in definition (fT6|) . while the theoretical construction in Section 
[3] used only the lowest lying states. 

5. Mapping to a thermalized system 

In the previous Section we defined the effective temperature starting out from the 
fluctuation-dissipation theorem, which is known to be valid for real thermodynamic 
systems at finite temperature. The idea behind our approach is a possible mapping of 
the subsystem in the infinite chain at T = to a finite system which is thermalized at 
some finite temperature. 

A more direct way of testing this mapping would be to compare expectation values 
in both of the above mentioned settings. That is, one has to consider an XX chain of 
finite length L connected to a heat bath, and try to adjust the temperature in such a 
way that the expectation value of a given quantity would be exactly the same as in the 
entangled block at zero temperature. 

Since the total transverse magnetization M of the finite chain is vanishing at any 
temperature for /i = 0, we shall choose the fluctuations to compare. Thus, the effective 
inverse temperature /3* is defined from the equality 

(M^)^. = (Ml), (33) 
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where the subscript on the bracket denotes an expectation value taken in the canonical 
ensemble with an inverse temperature /?*, while the right hand side is the T = 
fluctuation. We changed the notation of the effective inverse temperature to avoid 
confusion with former definition (fT6ll . The r.h.s of the equation was already determined 
in ( I20l) . thus one has to evaluate only the canonical expectation value. 

First we calculate the partition function Z = Tr(e~^^^). Since the XX Hamiltonian 
can be transformed into a set of noninteracting fermions, the partition function can be 
evaluated exactly for arbitrary values of the magnetic field: 

L 

Z{P,h) = l[{l + e'^'^), (34) 

k=l 

where = — cos(;^^) — h is the one-particle spectrum of H^. 

Starting out from the partition function one can calculate the moments of the 
magnetization as logarithmic derivatives respect to the magnetic field. The expectation 
value reads: 

^ fc=l ^ 

while the fiuctuations can be obtained as: 

(*^=>.-(M)?^^-^^-g^^^. (36) 

We have now a closed form for either side of the equation fl33l) . therefore (3* can 
be calculated numerically for arbitrary system sizes. Figure [6] shows the results of the 
calculations together with the asymptotic form obtained from the fiuctuation-dissipation 
ratio. Apparently, either definition for the effective temperature gives the same result, 
and the deviation (see inset) from the asymptotic curve tends to be relatively small even 
for moderate system sizes. 

The precise agreement of the inverse temperatures might look surprising, but it has 
a simple explanation in the given case. Namely it is the consequence of the temperature 
independent susceptibility. Indeed, looking at the formula of the magnetization ( |35l) . one 
has a sum of tangent hyperbolic terms, a function which is known to saturate for large 
arguments. For large system sizes f3* is expected to be large, thus the nontrivial terms 
in the sum are those where ~ 0, that is the Fermi level of the spectrum. However, as 
already shown in (fT2|) . the spectrum is linear near the Fermi level, thus application of 
a small field will cause a shift in the spectrum, nevertheless preserves linearity. Hence, 
the terms in the sum near the Fermi level cancel, the tangent hyperbolic being an odd 
function, and the result depends only on the difference in the number of terms added on 
the saturated tails. This is of course determined by the shift alone that is independent 
of temperature. 

Now, if the susceptibility is temperature independent and we use the definition (|T6ll 
of the effective temperature, the following relation holds: 

P*{M')p.=(3,ff{Ml), (37) 
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Figure 6. Effective inverse temperature obtained by solving Eq. (|33p numerically 
(squares). The solid line represents the asymptotic value Eq.([22|) obtained from the 
fluctuation-dissipation relation at zero temperature. Inset shows the relative deviation 
from the asymptotic curve. 



which, together with (l33l) results in /3* = Peff, 
identical. 



that is the different definitions are 



6. Subsystems in contact 



When one speaks about the temperature of a system there are some basic properties the 
definition has to fulfill. For example, when we take two systems with equal temperatures 
and place them in contact, the resulting system must have the same temperature as the 
initial ones. 

A possible analogue of this situation in our case is the following: we take two semi- 
infinite chains with blocks of length L on the boundaries and join them together. The 
resulting block of length 2L in the infinite chain has to have an effective temperature 
dictated by Eq. fl22l) . 

In order to solve this problem, we shall calculate the effective inverse temperature 
for a block of L spins at the boundary of a semi-infinite chain using again the fluctuation- 
dissipation approach for the transverse magnetization. Following the steps of Section 
14.11 we shall first calculate the magnetization with a nonzero magnetic field: 



(Ml) = Yl 



i=i 



TT 



9h 

dqsin^ (qj) - - 



TT Z 



l^sin(2jg.)^ (3^ 



TT 



J = l 



2j 
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where now brackets denote averages taken in the ground state of the semi-infinite chain. 

One can see that the first term of the above resuh coincides with the one obtained 
for the infinite chain ( fT71) . however the second term has an exphcit dependence on the 
location of the block sites which is the consequence of broken translational symmetry. 
Taking the derivative at zero magnetic field, the second term becomes zero for even and 
— l/vr for odd L thus for even L we recover the infinite chain result f|T8|) exactly. 

For the fluctuations we need the elements of the correlation matrix of the semi- 
inflnite chain which can be given as: 

sin(f(^-j)) sin(f(^ + j)) 

Now, one has to simply substitute Cij with C^^ in Eq. fl20l) . Introducing the difference 
r = i — j and the sum s = i + j of the sites we have: 

L 



2 L-l 2L- 



r=l s=2+r 
rodd sodd 



1 1 (-1)^ 

s rs 



4 («) 



We were able to evaluate the above sum only numerically. The result for the effective 
inverse temperature /3|jj(L) of a block of L spins in the semi- infinite chain is shown on 
Figure [7] together with the effective inverse temperature (3eff{2L) of a block of size 2L in 
an infinite chain. One can see an excellent agreement between the two results along with 
oscillations which might originate from the difference between even and odd block sizes. 
It is reasonable to think that these oscillations might be interpreted through the fact 
that the interaction between the spins generates anticorrelations between the transverse 
spin components, thus all quantities whose definition involves transverse correlations 
(such as the effective temperature) will be sensitive to the parity of the system size. 

In addition to verifying that the effective temperature does not change when 
bringing two subsystems of the same size into contact, we also obtained the formula 

which is valid when the block is at the boundary of a semi-infinite chain, and c is the 
constant defined in fl2T]) . This result also confirms the factors close to 27r instead of tt, 
obtained when fitting to the DMRG results with open boundary conditions. 

With the formula 0411) in hand, we can also confirm another thermodynamic 
property of the temperature. Namely, if we take two subsystems of different 
temperatures and bring them into contact the resulting temperature of the joint system 
must have a temperature between the two initial values. Changing to effective inverse 
temperatures, this means that the following inequalities should hold: 

/3:^/L)</3e//(L + L')</5:^;(L'), if L<L'. (42) 

In order to prove fj42l) . we use the equahty P^jj-^L) = Peffi^L) and require the positivity 
of the derivative of (3eff{L + L') respect to L and L'. This requirement yields the 
condition ln(L -|- L') + c > 1. As one can see, this condition holds even for the smallest 
possible unequal system sizes (L = 1, L' = 2) provided c > which is what we have 
found in all our calculations. 
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Figure 7. Effective inverse temperature of a block of 2L spins in an infinite chain 
(solid line) and a block of L spins on the boundary of a semi-infinite chain (circles). 
The inset shows the relative deviation from the /3e//(2L) line. The oscillations are due 
to the difference between even/odd block sizes. 



7. Final remarks 



We have presented an effective thermodynamic description of an entangled block in an 
XX spin chain at zero temperature. For a better understanding of the range of validity 
of this picture several comments are in order. 

First of all, it has to be noted that during our investigation of the one-particle 
spectrum of the effective Hamiltonian ([8]), we dealt only with the eigenvalues and did not 
look at the appropriate eigenstates. It was already demonstrated in [13] that, in contrast 
to the homogeneous plane wave eigenstates of the open XX chain, the eigenstates of the 
embedded block show enhanced amplitudes near the boundaries. Nevertheless, all those 
states with symmetrically occupied particle-hole excitations, such as the ground state, 
are identical. Furthermore, certain quantities such as the transverse magnetization have 
the same expectation values in either of the one-particle states. However, it is clear that 
the role of this strange amplitude enhancement requires deeper understanding, especially 
when dealing with quantities that are sensitive to boundary effects. It is possible that 
in such cases the validity of this simple thermodynamic picture might be lost, or at least 
requires more careful considerations. 

It would be also desirable to test the fluctuation-dissipation approach with some 
nonconserved quantity e.g. the longitudinal magnetization. Unfortunately, in such cases 
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the calculation of the response function is fairly involved, since the model becomes 
nonintegrable when the appropriate field is switched on. Even so, there remains 
the possibility to determine the response function numerically by means of DMRG 
calculations for finite chains. 

It is also reasonable to expect, that the effective temperature description might be 
carried over to other simple conformaly invariant models, such as the critical transverse 
Ising chain, where the form of the reduced density matrix of a subsystem is rather similar 
to the XX case [HI [lO], however the exact form of the effective temperature might be 
somewhat different. 

One can also speculate on the effect of placing the system in a heat bath of finite 
inverse temperature jShb- It is reasonable to expect that the effective temperature should 
not change dramatically as long as thermal fluctuations are subdominant of the quantum 
fluctuations, i.e. when Phb » Peff- On the other hand, when the temperature 
of the heat bath becomes much larger then the effective temperature << /3e//) 
the subsystem should be thermalized by the heat bath. Therefore, the behaviour of 
a subsystem should change, when the inverse temperature of the heat bath is on a 
scale (3hh oc This condition has been actually found for the transition temperature 
between the logarithmic and extensive behaviour of the entropy in gapless free fermionic 
systems [17], further supporting the existence of an effective temperature. 

As far as higher dimensional systems are considered, we remark on an other result 
of [T7] where the fluctuations of the transverse magnetization has been found to scale as 
L'^~^ InL with the system size in arbitrary dimensions. This property together with the 
extensive character of the susceptibility might lead us to speculate, that our fluctuation- 
dissipation approach would result an effective temperature with a same dependence on 
the linear size also in higher dimensions. We emphasize however, that these arguments 
require more careful investigations. 

Finally, it should be pointed out that there is an opportunity of extending our 
calculations to nonequilibrium situations, such as the driven XX chain. We believe, that 
the simple definitions we applied in the equilibrium case could be easily carried over, 
and the results may shed light on the possibility of introducing an effective temperature 
in systems far from equilibrium. 
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